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Abstract 

A free boundary problem arising from the optimal reinforcement of a membrane or from 
the reduction of traffic congestion is considered; it is of the form 

sup inf [ ( 1 ® | \7u\ 2 — fu] dx. 

f D 0 dx=m u ^ H o( D ) J D \ 2 J 

We prove the existence of an optimal reinforcement 6 and that it has some higher integrability 
properties. We also provide some numerical computations for 6 and u. 
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1 Introduction 


In the present paper we consider a variational problem related to various models arising in 
several fields among which shape optimization, optimal transport and elasto-plasticity. Here 
below we discuss the statement of the problem in some of them, the original motivation for the 
paper being the first one, for which we refer to and references therein. 

Optimal reinforcement of a membrane. Let D C R d be a bounded open set with a smooth 
boundary and let / G L 2 (D); D represents (when d — 2) a membrane, that we assume fixed 
at its boundary, and / represents a given load. The deformation u of the membrane is then 
obtained as the solution of the variational problem 


min 



/'«) 


dx : u G H, 


d (£>)}, 


where the elasticity coefficient in D is taken equal to 1, or equivalently of the PDE 


—A u — f in D, u = 0 on dD. 


We now want to rigidify the membrane by adding a reinforcement: this is modeled by a density 
function 6 > 0 which corresponds to changing the elasticity coefficient of D from 1 to 1 + 0. The 
total amount m of reinforcement is given and the ultimate goal is to increase as much as possible 
the rigidity of the membrane. This translates into the minimization of the elastic compliance, 
or equivalently into the maximization of the elastic energy, and we end up with the following 
optimization problem: 


sup inf [ ( ^ I \7u\ 2 — fu) dx. (1.1) 

f D 6 dx=m ueH^D)J D V 2 J 

Reduction to a traffic congestion problem. Consider D as an urban region with given traffic 
sources /; if a denotes the exiting traffic flux in D (or equivalently towards the boundary dD) 
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and if the cost of the traffic congestion in D is measured by the quantity f D ^\a\ 2 dx (see for 
instance m, 0 ), we may obtain a by solving the minimization problem 


. r [ 1, | 

mm < / -\<j\ 

^ Jd 2 


: dx 


div a = / in D j. 


Suppose that there is some given total amount m of resources that can be spent in order to 
decrease the traffic congestion: we assume that investing the amount 6 dx near a point xo 
produces a lower congestion, measured by the quantity ^(i+o) H 2 - Therefore, minimizing the 
total congestion cost leads to the problem 


mm 


1 


cr 


d 2 1 + 6 


dx : — diver = / in 


whose dual form is well-known to coincide with the optimization problem 0. 

The variational problem we deal with has been considered in the literature under various 
forms, below we mention two of them. 

The elastic-plastic torsion problem. Let D C K 2 be the section of a thin rod which is being 
twisted and let / be a constant corresponding to the twist of the rod. The function u G Hq(D) 
is a potential related to the elastic energy of the twisted rod. The response of some materials in 
this situation is of elastic-plastic type which means that the gradient of the potential u cannot 
grow beyond a certain threshold ft, i.e. u solves the variational problem 

mm {2 L \Vu\ 2 dx — f J udx : u G Hq(D), |V^| < ft|, 

the elastic region being given by {\Vu\ < ft} while the plastic one by {\Vu\ — ft}. 

In [2j Brezis proved that the potential u of the elastic-plastic torsion problem is in fact the 
solution of the problem 


— div ((1 + 9)Vuj — f in D, u — 0 on dD , 

where 6 is precisely the solution of the optimization problem 0 for some m > 0. In fact 9 in 
this case can be interpreted as the response of the material to the deformation of the rod. 

Several authors studied the regularity of the torsion function u and the free boundary 
5{|V^| < ft}. The optimal C 1,1 regularity of u was proved by Wiegner [Mj and Evans jH], 
while Caffarelli and Friedman provided a detailed study of the free boundary in [8j. 

Optimal transport problem In [Ij Bouchitte and Buttazzo considered a problem similar to 
0 without the underlying membrane, that is with 1+9 replaced by 0, and made the connection 
of the optimization problem 0 with the optimal transport theory. Several authors (see for 
instance m , m, d) studied the summability properties of the function 9 (called transport 
density ) in relation to the summability of the function / (the positive and negative parts / + and 
f~ are usually called marginals ). The usual setting considered in the literature is with D convex, 
Neumann conditions at dD , together with the zero average condition f D f dx = 0. In our case the 
Dirichlet boundary condition requires additional regularity of dD to get summability properties 
of 9. Moreover, the presence of the underlying membrane gives raise to a free boundary problem 
that we analyze in detail. 

The paper is organized as follows. The precise mathematical setting of the optimization 
problem and our main results are described in Section [2| while the proof is carried out in 
Section [3]. Section [5] is devoted to the elastic-plastic torsion problem , in which f — 1. Section 
[6] deals with some numerical computations providing the solutions u and 9 in some interesting 
situations. 
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2 Setting of the problem and main result 


Let D C be a given bounded open set and / G L 2 (D). For a given function u G Hq(D) and 
0 G L\(D) := G L l (D) : 9 > o}, the functional J f : H&(D) x L\(D) ^ M is defined as 

Jf(u, 6) = \ t f (1 + 9)\Vu\ 2 dx — ( uf dx, 

2 J D J D 


and the energy Ef{0 ) is given by 

E f (0) = min | J f (u,6) : u G F/q (£))}. (2.1) 

The minimum in 0 is always achieved due to the fact that every minimizing sequence is 
bounded in Hq(D) and so weakly compact in Hq(D). The minimizer is also unique due to the 
strict convexity of the functional Jf with respect to the first variable. If we denote it by uq we 
have that uq G Hq(D), \Vuq\ G L 2 (Q dx) and uq is the weak solution of 

- div ((1 + 0)Vue) = f in D, uq G Hq(D), (2.2) 

i.e. for every cp G Hq(D) such that |V(^| G L 2 (6dx) we have 


/ (1 + 6)Vuq • dx = / fipdx. 

Jd Jd 


Testing the equation (2.2) with (p = uq and integrating by parts leads to 

Ef(9 ) = Jf{u e , 9) = J^ugfdx. 


We consider the optimization problem 


max Ef(9 ), 

OeAm 1 


(2.3) 


where, for a given m > 0, the admissible set A m is given by 

A m = G L\{p) : J 9{x) dx < m|. 

Due to the monotonicity of Ef{9) with respect to 6 the formulation of the problem with the 
alternative admissible set 


A m = E L\(D) : J 9(x)dx = m^ 


is equivalent to the one with A m - 


Remark 2.1. It is interesting to notice that the optimization problem (2.3) is meaningful even if 
the load / is assumed to be a measure; in this case we have to allow 9 to be a measure too, and 
for every nonnegative measure 9 the functional Jf(u,9) and the energy Ef(9) are defined by 

J f (u,0) = \ [ (1 + 9)\Vu\ 2 dx - [ udf, VueCr{D) 

2 J D J D 

E f (9) = inf |j/(u,0) : ugC c °°(L»)}. 

If / is a measure of course we may have Ef{9) = —00 for some measures 9 ; this happens for 
instance when / concentrates on sets of dimension smaller than d — 1 and 9 is the Lebesgue 
measure. However, these “singular” measures 9 are ruled out from our analysis because of the 
maximization problem (2.3) we are dealing with. 
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Remark 2.2. Testing the energy Ef(Q) by the function u — 0 gives immediately the inequality 

E f (0) <0 V0 G Am. 

A more careful analysis (see for instance pQ) shows the better inequality 

i 2 (f) 


E f (0) < ~ 


2(\D\ + m) 


V6 € An 


where 


1(f) = sup [J udf : ueCr(D), |Vtt| < l}. 

On the other hand, if / G iL _1 (D) the energy Ef remains always bounded from below.: 

Jf(u,9)>- f (1 + 0)\S7u\ 2 dx - ||/||^-i(D)ll u lli?i(D) 

■' D 


l 


\^u\ 2 dx - U\\f\\ 2 H - 1(D ) + \\u\\ 2 h1o(d) ) > -i 


2 Jd ' 2 

Taking the minimum over u G Hq(D) we obtain 

Ef{9) > -x||/||#-i(D) 




V0 G v4 m . 


When / G L°°(D), setting M — ||/||oo and using the inequalities Ef(9 ) > Ef( 0) > £’m(0), 
which come from the maximum principle, we also have the inequality 


E f {0) > - 


/ wjjdx , 

Jd 


where wp is the solution of 


—A wd = 1 in D , wjj G (Z?). 


Our main result is the following. 


Theorem 2.3. Let D C 


6e a bounded open set satisfying the external ball condition and let 


f G L°°(D). Then there is a solution to the problem (2.3). Moreover, any solution 6 of (2.3) 
has the following properties: 

(i) (higher integrability) 6 G L P (D), for every p > 1; 

(ii) (min-max exchange) the energy Ef{6) satisfies 


1 


m , 


Ef(6) = min - / -||V^|| L oo — / uf dx; 

ueH^(D) 2 J D 2 Jd 


(2.4) 


(m) (regularity of the state function) t/ie solution uq of (2.2) is a/so t/ie minimizer of the 
right-hand side of (|2.4|). In particular, if f G C 2 (D) then uq G C fl,1 (D); 


(To) (optimality condition) 0 = 0 a.e. on the set {\Vuq\ < || Vi^Hloo}. 

In the case when the set D is convex we have a slightly stronger result. 

Theorem 2.4. Let D cK d be a bounded convex set, r G (d, +oo] and f G L r (D). Then all the 
conclusions of Theorem 2.3 hold. Moreover for the maximizer 6 we have that 6 G L r (D). 
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Before we pass to the proof of Theorem |2.3| and Theorem |2.4| we give some preliminary 
remarks in order to clarify the assumptions on D and / as well as to outline the main steps of 
the proof. 

The main difficulty in the existence result Theorem |2.3| is the lack of compactness in the space 
of admissible functions A m . Our approach consists in considering the optimization problem (2.3) 
on the admissible class 


An.p = {0 € L P (D), e>0, (J e p dxf /V < m}, 


(2.5) 


and then to pass to the limit as p —» 1. In fact, due to the weak compactness of the unit ball in 
L p (for p > 1) we can easily obtain the existence of an optimal reinforcement 6 p for 


ma x{Ef(9) : 9 E 

We discuss this problem in the beginning of Section |3| The rest of the proof is dedicated to the 
passage to the limit as p —>> 1. We notice that at this point we work with not just a maximising 
sequence of L 1 functions but with a family go maxima that obey some regularity properties we 
may use in order to get the desired convergence result. 

Our goal is to find a uniform estimate for 9 p in some Lebesgue space L r (D) with r > 1. In 
order to do that we use an estimate from m in a slightly more general form. Roughly speaking, 
it is of the form 


[ e r p dx<C r [ \f(x)\ r dx — [ F r (\Vu 0 \)Hqd 
JD JD JdD 

where C r is a constant depending on r, F r is a polynomial with positive coefficients depending 
only on r and Hqd is the mean curvature of dD. In the case when D is convex the boundary 
integral is positive, which immediately provides the uniform estimate which finally leads to the 
result of Theorem |2.4| On the other hand, if D is not convex then we should estimate the 


boundary term in a uniform way. At this point the bound from below of the mean curvature is 
natural to be imposed, but we have to estimate also the gradient of the solution on the boundary. 
We do this by constructing a barrier for uq p at each point of the boundary. In order to do this 
we use the external ball condition, which is stronger than just the bound on the mean curvature, 
and the assumption / E L°° which allows us to describe the barrier as an explicit solution in a 
suitably chosen domain. The construction of the barrier and the uniform estimate is proved in 
Section |3l 


In the case when f is a constant the proof of Theorem 2.3 can be obtained in a more direct 
way. In fact once the bound on the gradient is established on dD one can extend it in the interior 
of D by using the Payne-Philippin maximum principle HU- In this case also the boundedness of 
the optimal 6 is achieved. We notice that this result was already proved by Brezis in [2] where 
the optimal 6 is given through an explicit formula. 

Once we prove that the family 9 P is uniformly bounded in L r , we take, as a natural candidate 
for a solution, the weak limit of 9 P as p 1. In order to prove that the limit is a minimizer 
of (2.3), we use the characterisation of the torsion function u p as minima of functionals not 


depending only on / and p. 


There is an alternative approach to problem (2.3). In fact, the existence of an optimal 


reinforcement can be obtained for very general D and / if we relax the problem to the space of 
finite Borel measures on the compact set D. Moreover, for an optimal measure 9 the min-max 
exchange equality (2.4) still holds. Now under the assumption that D is convex one can use 


the regularity result for uq proved in [3] and then writing the equation for 9 as an optimal 
transport problem in which 9 corresponds to the density of transport rays one can obtain the 
absolute continuity and the summability of 9 using the results from m and m • We discuss 
this approach in the separate Section |4j 
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Remark 2.5. It is interesting to analyze the behavior of the solution u m of (2.4) asm4 +oo. 


mm 
veH^(D) 2 m 


Suppose that f > 0 and set v m = mu m . We have that v m solves the minimization problem 

^nJ D \ Vv \ 2dx + \\\ yv W 2 L™{D) ~ j D v f dx ■ 

The T-limit, as m +oo of the cost functional in the last line is 

^l|Vu ||J v fdx , 

hence v m = mu m converges in L 2 (D ) to the solution v of the minimization problem 

min hvvlliocm) - [ vj dx , 

VGH^(D) 2 y ’ Jd 


which is given by v(x) — Cr>ddD(x), where dQo{x) is the Euclidean distance from x G D to the 
boundary dD and Cd = f D doo(x) f (x) dx. 


3 Existence of optimal reinforcement 


In this section we carry out the proofs of Theorem 2.3 and Theorem 2.4 In order to make the 
presentation as clear as possible we divide the main steps of the proof into different subsections. 


3.1 Approximating problems: existence and regularity for p > 1 

In this section we consider the problem of finding an optimal reinforcement in the admissible 


class of reinforcements A m ,p given by (2.5). Precisely, for every p > 1 and / G L 2 (D), we 
consider the optimization problem 


ma x\Ef(9) : 9^Amv\ — max min Jf(u,9). 

OeAm, P ueHi(D) 

Interchanging the sup and the inf in ( |3.1| ) we obtain the inequality 

sup min Jf(u,0) < inf sup Jf(u,0). 

0^Am,p (D) u£zHq(D) 6^Am,p 


(3.1) 


(3.2) 


We now notice that in the Holder inequality f 9pdx < \\9\\lp the equality is achieved if 

and only if 9 = \\9\\lp Thus, for a fixed u G Hq(Q) the supremum on the right-hand 

side of (3.2) with respect 9 G A m:P can be computed explicitly: 


sup 

0^Am,p 


Jf(u,9) = - J |Vi/| 2 dx + ^ ^ J \S7u\ 2q dx^ — J fudx , 


where q—pj(p— 1) is the dual exponent of p > 1. If moreover u G W^ 2q (D), then the supremum 
is in fact a maximum and is achieved for 


9 = m \\/u\ 2 ^ q ~ l) ( J \Wu\ 2q dx ) 


-(«-i)/0 


Thus, for a given u , we can recover in a unique way the candidate for optimal 9 in (3.1) and thus 


it is now sufficient to optimize in u. We make this argument precise in the following proposition. 
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Proposition 3.1. For every p > 1 and f G L 2 (D) the optimization problem (3.1) admits a 
unique solution 9 p , given by 

d p = m\\7u p \ 2{ - q -^( K j \Vu p \ 2q d x y (q ~ 1)/9 , (3.3) 

where q = p/(p — 1) is the dual exponent of p and u p is the solution of the auxiliary problem 

minj^ J | Vu\ 2 dx + ^ ^ J \Vu\ 2q dx^j — J fudx : ueHq(D)^. (3.4) 


Moreover, if D C M. d is a bounded domain with C 1,a boundary and f G then u p G 

for a constant (3 depending on a, p, d, D and H/Hl 00 * In particular, 6 P is Holder 
continuous up to the boundary. 


Proof. As we already noticed, by exchanging the max and the min in the original problem (3.1) 
we get 

sup min Jf(u,9) < inf {- [ \X7u\ 2 dx+— ( [ \\7u\ 2q dx) — [ fudx\. (3.5) 

eeAm, P ueH'iD) neH^(D) 12 J D 2 \J D ) J D J f 

By the direct methods of the calculus of variations we have that the infimum in the right-hand 
side of (3.5) is achieved, the minimizer u p is unique and satisfies (weakly in W^ 2q (D)) the 


equation 


-Au p - mC p A 2 qU p = / in D, u p ^W^ 2q {D), 


where the constant C p is given by C p = ||^ and, for r > 1, A r is the r-Laplacian 
operator 

A r ^ = div(|V^r _2 Vi/). 


Defining 9 P as in (3.3), we get that u p satisfies 

-div ((1 + 9 p )Vu p ) = / in D, u p G Hq(D), 
in sense of distributions in D and so u p is the unique minimizer of Jf(-,9 P ) in L/q(D) and 

E f(9 P ) = Jf(u p ,9 p ) = ^ J \Vu p \ 2 dx^(J \X7u p \ 2q dx^j /q - j fu p dx , 

which together with (3.5) and the minimality of u p gives 

sup Ef{9) < Ef(0 p ). 

0 ^LsA-m ,p 

By the choice of 9 p we have f D 9 P dx = m p and by the fact that u p is the unique minimizer of 
(3.4), we get that 9 p is the unique solution of (3.1). 


The internal C 1,/5 -regularity of u p (and, as a consequence, the Holder-regularity of 9 p 


is 


classical and can be found in E51> while the up to the boundary version was proved in Theorem 
lof |Hj. □ 


3.2 The limit, as p —>► 1, of the state functions u p 

Let D C M. d be an open set of finite measure, / G L 2 (D ) be a given function, m > 0 and p > 1 
two fixed numbers; we denote by q = p/(p — 1) the dual exponent of p > 1. 

We define the functional F p : L 2 (D ) G R as 


F p {u) 


1 J \Vu\ 2 dx + y (y \Vu\ 2q dx} /q - J fudx, 

+oo, 


if ue W^ 2q (D), 
otherwise, 


(3.6) 
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while Fi : L 2 {D ) 


is given by 


F 1 (u)=\lI D ' Vu] 2 dX+ j' ]Vu]ll ~m-I D - 

[+oo, 


fudx , if G W’o’°°(-D), 

otherwise. 


(3.7) 


For p > 1 we denote by u p the unique minimizer of F p . 

Our main objective in this subsection is to show that u p converges to u\ in some suitable 
functional space (a natural candidate being L 2 (D )) and that the sequence of norms ||V^ p ||^ 2 g 
does not degenerate as p 1, a fact that we use later in the uniform estimate of ||0 p ||^oo. Before 
we continue with the result from this section we recall the notion of T-convergence introduced 
by De Giorgi: 

Definition 3.2. Let X be a given metric space. We say that a sequence of functionals F n : 
X — y 1R. U {+oo} r -converges to the functional F : X —>► R U {+oo} ; if 

• (r — lim sup) For every sequence x n G X converging to x G X we have 

F(x) < lim inf F n (x n ); 

• (r — lim inf) For every x G X there is a sequence x n G X converging to x and such that 

F(x) > lim sup F n {x n ). 

Proposition 3.3. Let D C R. d be an open set of finite measure and let f G L 2 (D). With the 
notation introduced in (3.6) and fl3?fl ) we have: 

(i) The sequence of functionals (F p ) p> i T-converges in L 2 (D) to the functional F\. 

(ii) The sequence of minima {u p ) p> i converges strongly in H l (D) to u\. 

(Hi) lim ||V^J| L 2 q = \\Xui\\l°°- 

p—> l 


Proof. The proof is well-known, let us give it for the sake of completeness. Suppose that 
v p — y v in L 2 (D ) and lim F p (v p ) < +oo. 

p—>i 

Then there is a constant C > 0 such that for p > 1 small enough we have 

^||Vv p ||| 2 + - J^Vpfdx < C. (3.8) 

Thus, by the Gagliardo-Nirenberg-Sobolev inequality we have that v p is uniformly bounded in 
Hq(D) and so we can suppose that v p converges weakly in Hq(D) to v G Hq(D) (which gives 
the semicontinuity of the first term of (|3.8|)). We now notice that 


\\\7v\\ L oo(d) = sup| J Xv'vdx : v G C(P(D; R d ), J \v\dx<l^. 

On the other hand for a fixed v G C^°(D; R rf ) with f D \ v\ dx < 1, we have 


/ Xv • is dx = lim 

Id p^°°Jd 


< 


Xv p • v dx 


( I \vv p \ 2 u x y /(2q \j in w*- 1 ) dx 

< (liminf ||Vu p || L 2,^ ( lim |M| L 2 P /(„+i) j 


(2q-l)/(2q) 


< \\n\\ L i lim inf II 11 2,2 9 , 

p->i 



which proves the semicontinuity of the second term in (3.8). 

In order or prove the T — limsup inequality consider a function v G and take v p to 

be the constant family v p = v, for every p > 1. Since in this case we have || \7v\\ L 2 q || 

as q —» +oo, we obtain 

F(v) = lim F p (v), 

p—>i 


which concludes the proof of (i). 

We now prove (ii). Since u p is the minimizer of F p and F p ( 0) = 0, we have that F p (u p ) < 0, 
which gives 


\\\^ u pW 2 l^(d) - J D u pf dx < II u pIU 2 II/IIl 2 < M D ) 1/2 IIV« P |U2||/|| L 2, 

which gives the uniform (in p) bound || \7u p \\ L 2 < 2 \\f\\ L 2 \i(D)~ 1 G, Thus, there is a subsequence 
of u p converging in L 2 to some u G L 2 (D). Since F p T-converges to Fi, u is necessarily the 
minimizer u\ of F\ (and since for every subsequence the limit is the same u p converges to u\ in 
L 2 ). Moreover, by the optimality of u p and the T — lim inf inequality we have 

hmsupFp(?ip) < liminf F p {u\) = F\(ui) < liminf F p {u p ), 
which together with the fact that f D u p f dx —» f D uf dx gives that 

||Vui ||£ 2 +m||Vwi|||oo = lim |'||V«p ||^ 2 +m||Vitp||| 2 ,|. 

Now since both the terms are semicontinuous with respect to the L 2 convergence of u p to u\ we 
get that 

||Viti||| 2 = lim ||Vu p ||| 2 and ||Vni|||oo = lim ||Vu p ||| 29 . 

p—>l p—>l 

The first equality proves that u p u\ strongly in H l (D), while the second equality is precisely 
(in). □ 

3.3 Construction of a barrier for u p 

In order to obtain uniform estimates for u p and 6 p we impose some mild geometric assumptions 
on the domain D. Precisely we assume that D satisfies a uniform external ball condition that 
is there is some p > 0 such that for every point of the boundary xo G dD there is a ball 
B p (yo) C R. d \D such that xo G dB p (yo). This assumption allows us to construct a barrier for u p 
and thus to give an estimate on the gradient | \7u p \ on the boundary dD in terms of the diameter 
of D and the largest possible p. In the case when D is convex we have an estimate only in terms 
of the diameter diam (D). In this estimate we use a weak maximum (or comparison) principle 
which we state below in a general form. 

Lemma 3.4 (Weak maximum principle). Let D C R. d be a 
[0,+oo) —>► [0, Too) be a convex function such that G f ( 0) > 0. 
denote by uqj the unique minimizer of the problem 

min|- f G(| V?i| 2 ) dx — f uf dx : u G 

^2 J D Jd 

(a) If f,g G L 2 (D) are such that f > g, then unj > 

(b) If uo dD is a bounded open set and f > 0 a.e. on D, then uqj > u u j. 

We are now in a position to construct our barrier. 


bounded open set and let G : 
For a function f G L 2 (D) we 
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Lemma 3.5 (Pointwise boundary estimate). Let D C R. d be a bounded open set whose boundary 
dD is locally a graph of a C 1,a function for some a G (0,1). Let xq G dD and let B p (y o) be a 
ball in R. d \D tangent to dD in xq. Then, the minimizer u p of F p satisfies the gradient estimate 


(l + mCp\S/up(xo)\ 2 ^ q ~ 1 ^ |V« p (a;o)| 


< 11 / \ \L°°(D) 


1 diam(i/) \ d 1 


P ) 


diam (D). 


(3.9) 


Proof. Without loss of generality we assume yo — 0. Let R be the diameter of D. Then 
D C Br+p \ B p . We set for simplicity 


M — 


loo , C — mCp — m(^J \S7u p \ 2q dx'j 


, -!/p 


and we stress that M and C are fixed constants in the proof below. Let u be the solution of the 
equation 

-div ((1 + C|Vu| 2(<z_1) )V«) = M in D, u G Wq’ 2q (D), 
and let U be the solution of 

-div((l + C|Vf/| 2 ('?- 1 ))VC/) =M in B r+p \B p , U € W^ 2q (B R+p \ B p ). (3.10) 

By the maximum principle, applied to the operator associated to the function G{t) — tpCt q ~ x , 
we have \u p \ < u < U. On the other hand, the solution U of (3.10) is radially symmetric, i.e. a 
function of r that we still denote by U and that solves 

-r l - d d r (r d - l (l + C\d r U\ 2 ( q -V)d r lP) = M in }p,p + R[, U(p) = U(p + R) = 0. (3.11) 


Let pi G]p, p + R[ be a point where U achieves its maximum. Integrating (3.11) we have 
p d - 1 {l + C\d r U{p)\ 2{q - l) )d r U{p) = [ P1 Mr d ~ l dr < MR(R + p ) d ~ l , 


which gives (3.9) since d r U(p ) > 0. 


□ 


Taking the largest possible p in the pointwise estimate (3.9) we obtain the following result. 

Lemma 3.6 (Boundary estimate for a regular domain). Let D C M. d be a bounded open set 
whose boundary dD is locally a graph of a C 1,a function for some a G (0,1). 

(i) If D is convex, then we have 


1 + mC p \\S7u p fiif D) ) ||Vn p || L oo ((9D) < diam(L»)||/|| L oo (£>) . 


(ii) If D satisfies the external ball condition with radius p, then 

(l + mC p \\Vu p \\l { rf D) ) \\Vu p \\ Loo{dD) < ll/llicop) diam(D) (l + . 

3.4 Uniform estimate for 9 P 

In this section we extend the boundary estimate for 9 P to the entire domain D. We need 
a uniform estimate on the norm ||%,||loo, which is classical and so we give only the precise 
statement here. 
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Lemma 3.7. Let f G L r (D ) for some r > d/2. Then we have the estimate 

where Cd is a dimensional constant. 


We now give our main a priori estimate for 8 p . The statement might appear as a generaliza¬ 
tion of an estimate by De Pascale-Evans-Pratelli, but the proof is precisely the one they gave in 
m We reproduce the proof below for the sake of completeness and in order to show that the 
presence of a general functional G and the non-convexity of the set D do not influence the final 
result. In fact the only difference is that since the domain is not convex there is a boundary 
term that appears in the final estimate. 


Lemma 3.8 (De Pascale-Evans-Pratelli). Let D be a bounded open set with smooth boundary 
and let f G L r (D ) for r >2 and r > d/2. Let G : [0, +oo) —>> [0, +oo) be a convex function such 
that G'( 0) > 0 and let u G C 1 (D) H H 2 oc (D ) be the solution of 

— div (G'(\ \7u\ 2 )\7u) = / in D, u = 0 ondD. (3.12) 


Then for every £ G (0,1) we have the estimate 

[ \G'(\Vu\ 2 )\ r \Vu\ 2 dx<3£ [ \G\\X7u\ 2 )\ r dx 

JD JD 


+ 


(r — l) r 
c-r—l 


+ 


1 


,2r-l 


\U\\ 


r 

L r 


(r-l) 2 |M|| 0 


£ JdD 

where H is the mean curvature of dD with respect to the outer normal. 


[ H\G\\Vu\ 2 )\ r \Vu\ 2 d'H d -\ 

JdD 


Proof. We first notice that under the assumptions that D and / are smooth, we have that u is 
also smooth up to the boundary dD. Thus we can use the multiplication technique from m in 
order to obtain the L p estimate of G'dV^I 2 ). For the sake of simplicity we set a — G'dV^I 2 ). 
We first test the equation Q3.12| ) with the function a r ~ 1 u G Hq(D) obtaining 


Id 


a r \Vu\ 2 dx + (r — 1) f ua r • Vcr dx = f fua r 1 dx 
J D Jd 


< IMIL° 


L 


L r \ I V r dx 


1 — 1/r 


(3.13) 


Our main objective is to estimate the second term of the left-hand side. Multiplying (3.12) by 
div(cr r_1 V^) = ( a r ~ l Uj)j we get 

f (aUi)i(a r ~ 1 Uj)j dx = — f f(a r ~ 1 Uj)jdx = — f f(cr r ~ 2 auj)jdx 

Jd Jd Jd 

— — f fa r ~ 2 (aui)i dx — (r — 2) f fcr r ~ 2 <jjUj dx (3.14) 
Jd Jd 

< [ f 2 a r ~ 2 dx + (r — 2) [ \f\ a r ~ 2 \ Va • \7u\ dx. 

Jd Jd 

On the other hand, we can integrate by parts the left-hand side getting 
/ (aui)i(a r ~ 1 uj)j dx = — / aui(a r ~ 1 uj)ji dx+ ( o- r ~ l Uj)j(jUini dH° 

Jd Jd Jod 


jd -1 


= / ( aui)j(a r 1 Uj)idx-\- / (cr r 1 Uj)jaUiUi — aui(a r 

Jd JdD L 


dU 


d-i 


= / ( aui)j(a r 1 Uj)idx-\- / a r (ujjUiUi — UijUiUj) dTL 


d— 1 


id 


IdD 


/ (aui)j(a r 1 Uj)i dx + / a r u n (Au - u nn ) dH d 1 , 
Jd Jod 
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where u n and u nn indicate the first and the second derivatives in the direction of the exterior 
normal n to dD. Developing the term in the volume integral we have 



/„(- 


'\\Hess(u)\\i + (r - l)a r ~ 2 \Va • Vu\ 


+ ra r 1 ajUiUij^ 


dx 


+ [ cr r u n (Au — u nn ) dTL d 1 , 
J dD 


where \\Hess(u )\\2 = Ylij u ij — 0* We first notice that aj = 2ukUkjG"(\ V^| 2 ) and thus, by the 
convexity of G we have ajUiUij > 0. On the other hand, on the boundary dD we have 


Au — U nn “1“ H-Ujn 


where H is the mean curvature. Thus we have the inequality 

[ (o‘Ui)i((j r ~ 1 u j ) j dx > [ (r — l)cr r-2 |Vcr • X7u\ 2 dx + [ Ha r \Vu\ 2 dU d ~ x , 
J D J D JdD 


which together with (3.14) gives 

[ a r_2 |Va • \7u\ 2 dx < [ \f\ 2 a r ~ 2 dx- [ Ha r \X7u\ 2 dU d ~ x 
Jd Jd J dD 

< WfWlAWt- 2 ~ [ Ha r \S7u\ 2 dH d - 1 . 

JdD 

Using (3.13) we can now repeatedly use the Young inequality A a B@ < zolA + z a ^/3B for 
0 < cq (3 < 1 and a + (3 = 1 obtaining 

f (j r \\7u\ 2 dx < (r — 1) f (j r—1 H^Hl 00 | Vi/ • Vcr| dx + H^Hl 00 \\f\\L r ll^lli^ 1 
Jd Jd 


<e [ a r dx+ ^ ^ [ a r ~ 2 \Wu • Va | 2 dx + ||u|| 

Jd £ Jd 


II _ IIr —1 

L r 11 & 11 J^r 


< £ / cr r dx + 

Jd 


(r — l) 2 ||id |2 


i-IMir 2 + IMMI/IUHMlU 


(r — l) 2 ||?i 112 


I L°° 


[ H (j r | V 1 
JdD 


u\ 2 dU d ~ l 


<Zedx+ (U^zr- + ^=r)l 


(r — 1 ) 2 || u 112 


\L°° 


[ Ha r |V 

JdD 


U Wl°°\\J II Lr 

u\ 2 dTL d ~ 1 , 


which is precisely the claim. 


□ 


In order to give a uniform estimate for 9 p on every domain satisfying the external ball 
condition and not only on the regular ones, we use an approximation argument. The following 
lemma that we use is well known in the 7 -convergence theory. 

Lemma 3.9. Suppose that D cK d is a bounded open satisfying the external ball condition and 
that D n is a sequence of open sets of finite measure containing D and such that \D n \ D\ 0. 
Let f G L 2 (Di), p G [ 1 ,+ oo ) be fixed and let u p be the minimizer of F p on D and u p on D n . 
Then the sequence u p converges to u p strongly in i7 1 (R d ) and W 1,2q (R d ). 

Proof. By the Sobolev inequality it is immediate to check that the sequence vA n is uniformly 
bounded in H 1 (and lD 1 , 2 <? (M d ) and so we may suppose that up to a subsequence it converges 
in L 2 to a function u. Moreover, all the functions u p are uniformly bounded in L°° and so, the 


12 










function u is zero almost everywhere on R d \D. Now, by [?] we have that u G Hq(D). Moreover, 
for any function v G Hq(D) C Hq(D u ) we have 

i^(V) < liminf F p {u p ) < F p (v), 

which proves that u = u p is the minimizer of on D. Moreover, since F p (u) > F p (u p ) for every 
n, we get that F p {u p ) -+ F p (u). By the strong L 2 convergence of u p to u p we obtain also that 
||V ^||^2 -+ HViXpll^ and ||V^||^ 2 g -+ ||V^||^ 2 g , which concludes the proof. □ 

Proposition 3.10. Let D C R rf be a bounded open set of finite perimeter satisfying the external 
ball condition and let f G L°°(R d ). Then, for every r > d, there are constants 5 > 0, depending 
on D and f, and C, depending on r, on the dimension d, the perimeter P(D), the diameter 
diam(D), the radius R of the external ball, the norms ||/||l oo ( j d) an d H^iHl 00 ; such that 

\\O p \\ L r(D) < C for every pG (1,1 + 5). 

Proof. Suppose first that D has smooth boundary. We recall the notations u p for the minimizer 
of the functional F p in Hq(D), 6 p for the optimal reinforcement, and C p — ||Vi^|| 2g 2 ^ where 
Q — p/(p— !)• Setting 

G p (i) = * + !++ 

we have that u p is the minimizer of the functional 

Hl(D)3u^~ [ G p (\S7u\ 2 ) dx - [ ufdx. 

2 J D J D 


Moreover, we have 


0 p = m( ) =mCp\Vu p \ 2( - q ^ = G' (\Vu p \ 2 )-1. 


By Lemma 3.8 and the mean curvature estimate Hqd > — l/i2, we have that 

[ cr r p \\7u p \ 2 dx < 3e [ dx + +-+-) 

jd Jd \ e e j 

(r-i) 2 Kii|« r 
eR J a 

We denote by S and B the sets 


JdD 


v\\l° 


cr r p \\7u p \ 2 dU d ~ X , 


(3.15) 


S = {x G D : \\7u p (x)\ < \\\7u p \\ L 2 q } and B = {x G D : \Vu p (x)\ > \\Vu p \\ L 2 q }, 


and we notice that we have the inequality a p < 1 + m on S, where we set for simplicity 
a p — G p (\\7u p \ 2 ). We now estimate separately the three terms on the right-hand side. For the 
first one, after integrating separately on S and B , we get 



D dx+ L 


a. 


\dx < (1 + m) r \D\ + ||Vu, 


\-2 
I L 2 « 


L 


a n 


|| \7u p \ 2 dx. 


(3.16) 
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We can estimate the second term simply by using Lemma [3T7| while for the third one we have 

[ a r \Wu p \ 2 dU d ~ x = [ <7p|V-u p | 2 d'H d ~ 1 + [ a r \Vu p \ 2 dTi^ 1 

JdD JSndD JBndD 


< (1 + m) 


r —2 


/ SHOD 


a p \Vu p \ 2 d'H 


d— 1 


jd -1 


+ 11 / a r p \Wu p \ r dn° 

JBndD 

< (1 + m) r ~ 2 P(D)\\f \\loo diam (D) 2 (l + dia ™( jP ) 

diam(P)\ r ^ _1 ) 
P ) 


+ ||V^||^;P(D)||/|r L oo diam(P) r 1 + 


Now taking in (3.15) e — W'V u p \\ 2 L2q / 6 and absorbing part of the first term in the right-hand side 
of ( 3.15| ) into the left-hand side we get 

\ f D a P^ Up \ 2 dx - ^H^plli 2 ^ 1 +™) r \ D \ 

(r — l) r ( 1 \ 

JnF 2 ) 


+ 6 


2r—l 


(i 


V Unr 


\2r—2 


+ 


11 L° 


+ 


6 pIIl 2 ^ 
6 (r - l) 2 ||n p |||oo 

II 


|| VMpll L 2 9 

(l + m) r - 2 P(D)\\f\\ 2 LOO diam(£>) 


2 ^ + diam(L>) ^ 2(d 1} 


+ 


iivu p ii^p(d)h/ii£„ diam(D)-' y + di y p > y (J 11 


We now notice that the same inequality holds for irregular D. In fact we can approximate 
D by a sequence of smooth sets D n containing D such that | D n \ D\ -G 0, P(D n ) < 2P(D), 
diam(P n ) < 2diam(P). By Lemma 3.9 the right-hand side of the above estimate passes to the 
limit, while the left-hand side is semicontinuous which finally gives the estimate passes to the 
limit. We now consider S > 0 such that 

1 3 

dl^ilk^CD) < IMIl°°(.d) < ^\\ui\\l^{d), for every p<E(l,l + < 5 ). 

Then there is a constant C = C(d,r,m, ||/||l°°, P(D), diam(D), ||'i*i Hz, 00 ) such that 


Jd 


<j r p\\7u p \ 2 dx < C , for every p E (1,1 + S). 


Now the conclusion follows by the inequality f D 9 p dx < f D a r p dx and by applying one more time 
the estimate (3.16). □ 


Uniform L°° estimate for 9 p in the case of constant force. In the case / = const on D 
a uniform estimate on 9 p can be obtained in a more direct way by using the Payne-Philippin 
maximum principle m that we recall below. The idea is similar to the one used by Kawohl in 
m for the infinity Laplacian. 

Theorem 3.11 (Payne-Philippin maximum principle |18|). Let D cR rf be a bounded open set 
with C 2+£ boundary and let g : [0, +oo) -G [0, +oo) be a C 2 function such that 2 sg'(s) + g(s) > 0 
on [0,+oo). Then, the function 

r\S/u{x)\ 2 

G(x ) = / (2 sg'(s) + g(s )) ds, 

Jo 

achieves its maximum on the boundary dD, where the function u is the solution to the equation 
— div (p(| \/u\ 2 )Vu) = 1 in D, u — 0 on dD. 
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We are now in a position to give a uniform estimate on ||^||l oo ( j d)- 

Proposition 3.12. Let D C R. d be a bounded open set satisfying an external ball condition with 
radius R G (0, +oo], where by R — +oo is intended that D is convex. Then, there is e > 0 such 
that for p G (1,1 + e] we have the estimate 

Op < 2||Vui||“ 2 (D) M 2 diam(T>) 2 (1 + dia ™(-°) ^ 


For convex D we have 

Op < 2||V«i||^ 2 M 2 diam(D) 2 . (3.17) 

Proof. We consider only the case of D convex since the two cases are analogous. We first prove 
(3.17) for regular domains and we then argue by approximation. Indeed, suppose first that D 
has regular C 2,6 boundary. Then u p is a solution of 


— div {g p (\Vup\ 2 )'Vu p ) — M in D, u p = 0 on dD , 


where g p is given by 


g p (s) = s + mC p s q \ 


and is sufficiently smooth when q > 1 is large (i.e. p close to 1). By the Payne-Philippin theorem 
we have that the function 

G p {x) = h\7up\ 2 + — — -mCp\Vu p \ 2q , 

Z q 


assumes its maximum on the boundary of D. On the other hand, by (3.17) on dD we have 
2q-l. 


G p < 


-\\7u\ 2 (l + mCp\\/u\ 2( - q ~ l) ^ < — —^M 2 diam(H) 2 , 


Q 


and so on the entire D we have the estimate 


mC p \Vu\ 2q < ^ G p < M 2 dmm(D) 2 , 

ZjCf _L 


The optimal potential 0 P satisfies 


01 


= O^i-V = (mCp\\7u p \ 2 ( q -Vy iq 1} = (m^Cy^-^mCplVu, 


1 2 q 


We now recall that by the definition of C p we get 


cyb-v = (wvuj-y p ) 1/iq 1} = ii vu "- 2 


T>N.L 2 9 ) 


Thus we obtain 


0 P < (ro p/<? ||Vu p ||“ 2 2 M 2 diam(T>) 2 ) 


'P\\L 2c r 


1/p 


The case of irregular D follows by approximation with smooth convex sets containing D and 
then applying Lemma |3.9| □ 
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3.5 Conclusion of the proof of Theorem |2.3| and Theorem |2.4 

Proof of Theorem \2.3\ Let for p > 1, u p G W$ q (D) be the minimizer of F p and let 9 p be the 
optimal potential relative to u p . By the fact that 6 p G A m , P and by Proposition 3.10 we have 
that there is a constant C such that for p > 1 small enough we have 


Id 


9 P dx — m p 


and 


L 


e r v dx < c, 


for some r > d. In particular, 9 p is uniformly bounded in L 2 (D) and so, up to a subsequence, 
9 P converges weakly in L 2 (D) to a nonnegative function 9 G L 2 (D). Moreover, we have 


[ 9 dx = lim [ 9 P dx < liminf ||0 j9 ||^p| J D| 1 / g = 

JD p-^-l Jd P-+1 


m. 


On the other hand, by Proposition 3.3 we have that as p —» 1 the solutions u p converge strongly 
in Hq(D) to the minimum u\ of the functional F\. Thus, for every fixed smooth function with 
compact support p G C%°(D) we have 


/ fpdx— / (1 + 9 p )\7u p • \7p dx - 

Id Jd p- >1 


[ (l + 0)Vui 

JD 


Vp dx. 


Thus, we have that u\ is in fact the solution of the equation 

— div ((1 + 0)Vuij — f in D, u\ = 0 on dD , 
and after integration by parts 

E f(0) = \ f (1 + 0)\Vui\ 2 dx - f fuidx = -\ [ fuidx. 
z Jd Jd z Jd 

By the convergence of u p to u\ in L 2 (D) and by Proposition 


3.3 


we have 


Ef(6) = ~- 
which proves that 


[ fu\ dx — lim — [ fu p dx = lim Ef(9 p ) = lim F p {u p ) = F\(u{), 
Id p ^ 1 2 Jd p v 


Ef(9 ) = min FUu). 

ueH'(D) 


(3.18) 


On the other hand we have the general min-max inequality 

sup Ef{9) = sup min Jf(9,u) < min sup Jf(9,u) = min F\{u), 

9^A.m O^Am u £Ho(D) OElAtti uEHq (D) 

which concludes the proof that 9 is a solution to the problem 0- 

Claim (i) follows by the fact that 9 G L r (D ) and that r can be chosen arbitrary large. 
Claim (ii) is in fact (3.18) and the claim (iv) follows directly by this equality. 

The regularity of the minimizer u\ of F\ (in) was proved by Evans in [13]. □ 

Proof of Theorem \2A\ We first consider the case r < +oo. We notice that in the case of a 
convex set there is no boundary term in the De Pascale-Evans-Pratelli m estimate. Thus, for 
/ G L r (D ) there is a constant C r such that for every p large enough 

11 11 i/ r < C r (l + ||/||i/ r ) • 

The rest of the proof follows as in the proof of Theorem |2.3| □ 


Remark 3.13. By using methods of optimal transportation theory slightly finer results can be 
obtained m , private communication): 

- if D has the exterior ball condition and / G L P (D ) with p G [1, +oo], then 9 G L P (D)] 

- if D is any bounded open set and / G L p oc fD) with p G [1, +oo], then 9 G L p oc (D). 
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4 Alternative approach to the optimal reinforcement problem 


Let D C be a bounded open set with a smooth boundary. The right hand-side / represents 
a given force on the membrane D. We allow / to be a general signed measure / G Ai(D). A 
reinforcement of the membrane is given by a measure fi in the class of nonnegative 

measures on D. For a vertical displacement u G C^(D), the energy of the membrane subjected 
to the force / and reinforced by (i is given by the functional 

Jf(u,pL):=- f \Vu\ 2 dx-\ — f {Vu^d/i — f udf , 

2 Jd 2 J D J D 

and we define the Dirichlet energy functional 


Ef(n) := inf j Jf(u, pf) : u E C^{D) j. 


As we noticed in Remark 2.1 for some measures / and pi we may have Ef(pi ) = — oo. This occurs 
for instance if the force / ^ H~ l {D ) and (i is the Lebesgue measure on D. However these cases 
are ruled out from our discussion, because we are interested in the maximal reinforcement of the 
membrane. In fact, for any measure / there exists a measure pi, corresponding to the transport 
density of a suitable Monge optimal transport problem, such that Ey(/i) > — oo (see for instance 
[1]). The optimization problem we are interested is the following: 


max 


: [i G *4.77^, 


(4.1) 


where the parameter m > 0 represents the available quantity of reinforcement material and the 
admissible set Am is given by 

A m = |/x <E M + (D) : n(D) < mj. 

A first existence result in the class of nonnegative measures fi is the following. 

Proposition 4.1. Let D C M. d be a bounded smooth open set and let m > 0. Then for every 


signed measure f G M(D) the optimization problem (4.1) admits a solution (i G A i + (D) with 
jd(D) = m and jd(dD) = 0. Moreover, we have that 


E f (p) = in£ _ 

u£Wq'°°(D) J 


\L |vt 


m , 


' dx + —\\Vu\\ L °o 


-L 


udf. 


Proof. We first notice that for every fixed u G C^{D) the map /x i-G Jf(u,/i ) is continuous for 
the weak* convergence. Hence the map pi i—>> Ef(/i ), being an infimum of continuous maps, is 
weak* upper semicontinuous. The conclusion follows by the weak* compactness of the A m and 
the fact seen above that there is at least pi such that Ef(pi) is finite. The last claim follows 
by the monotonicity of the functional Ef and the fact that for every measure /a G A m it holds 

To prove the last claim we notice that the functional Jf : A m x Cl(D) —>> M satisfies the 
assumptions: 

• Jf is concave and upper semicontinuous in the first variable with respect to the weak-* 
convergence; 

• Jf is convex in the second variable. 

According to a classical result (see for instance m we may interchange the inf and the sup for 
Jf obtaining 

sup inf Jf(pb,u)= inf sup (4.2) 


neAm u ^Cl(D) 


ueCl(D) ^eArr 
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The supremum with respect to fi is easy to compute and we have 


sup = [ |Vix| 2 dx + ^-||V^|||oo — f udf 

/^iG>A.rn " JD " JD 


as required. □ 

Therefore, we end up with a variational problem for the functional 

Fi(u) = \j D | Vn| 2 dx + j\\Vu\\l -J udf. (4.3) 

The existence and the regularity of the minimizers of F± has been widely studied in the literature 
(see for instance [1] and HU)- We summarize the known results below. 

Theorem 4.2. For every measure f G Ai(D) the optimization problem 

min {u(0 = ueH^D)} 
admits a unique solution u. Moreover, 

(i) if f G L P (D) with 1 < p < +oo ; then u G W 2,P (D ); in particular if p > d, then the solution 
u is C 1,a (D) for some a > 0; 

(ii) If D is convex and f G L°°(D ) then A u G L°°(D ); 

(in) if f G C 2 (D), then u G C l ^{D). 

Proof. The existence follows by the direct methods of the calculus of variations while the unique¬ 
ness is a consequence of the strict convexity of the functional F±. The solution u is clearly 
Lipschitz and setting M — HV^Hl 00 we have that u is also the solution of the problem 

minj- J \\7u\ 2 dx — J udf : u G \\\7u\\l°° < M j. 

The claims (i) and (ii) were proved in [2], while (Hi) was proved in [2]. □ 

Let now / G L P (D) with p > d and let jl be the solution to the problem ( |4.1| ). We notice 
that •) can be extended from a functional on C^(D) to a functional on (m[D ) n Cq(D). 

Thus, we can use the minimizer u as a test function in Jf{p, •) obtaining 

F\{u) > Jf%u) > inf _ Jf(fi,u) = E f (fi ) = Fi(u), 

ueC 1 (D)nC 0 (D) 

where the last equality follows by ( |4.2| ). Thus the second inequality is an equality and implies 
that u is the minimizer of the functional •) in C l (D) D Co(D) and so it satisfies the Euler- 
Lagrange equation 


— div ((1 + fxjVu) — f in D, u = 0 on dD. 
Using u as a test function in ( |4.4| ) we get 

f \S7u\ 2 dx-\- f \\7u\ 2 djl— f uf dx. 

JD JD JD 


(4.4) 


On the other hand, since u is a minimizer of F\, we have that the function t F\(tu) has a 
minimum in t — 1, which gives 



dx + m\\X7 u\\ 2 L oo = 


L 


uf dx. 
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Putting together the two identities we have that 


I \Vu\ 2 dp — m\\S7u\\ 2 L oo — mM 2 , 
Jd 


which together with (|4.4|) gives the Monge-Kantorovich transport equation 


— di v(/iV^) = / + A u 
\Vu\ < M, 

|V4Z| = M on spt(/i). 


in D, u = 0 on dl?, 


(4.5) 


This equation was widely studied in the framework of optimal transport theory. In particular, 
it was proved in m, Em and m that the integrability properties of the right-hand side give 


the integrability of p. Thus we can prove directly Theorem |2.4| using these optimal transport 
results. 

Alternative Proof of Theorem 


2.4[ Let / G L r (D ) for some r > d. By Proposition |4.1| there is 
a solution ft to the problem ( 4.1| ). By the regularity result of Brezis and Stampacchia [3] the 
minimizer u of F\ is C 1,a for some a > 0 and A u G L r (D). Thus jd and u solve (4.5), whose 
right-hand side is in L r (D). Now by the summability results from [12] and m we have that 
li G L r (D ), which concludes the proof in the case r < oo. If / G L°°, then again by [1] we have 
that A u + / G L °°. Now by jT9] (and also [20]) we have that also jd G L°°(D). □ 


5 The elastic-plastic torsion problem 

The solution of the optimization problem 0 is related to the solution of the variational 
minimization problem (4.3). In fact, the solution of (4.3) is precisely the torsion function u 


associated to the optimal elasticity term 6. On the other hand, in the case / = 1, the problem 
(4.3) is equivalent to the well-known of elastic-plastic torsion problem. In fact if u G VF 1,00 (D) 


is the solution of the variational problem (4.3), then setting n — HV^Hl 00 , it is immediate to 
check that u is also the solution of 

minj^ J \\7u\ 2 dx — J udx : u G Hq(D), ||V?/||lo° < ^j. (5.1) 


This problem was a subject of an intense study in the past. The first results are due to Ting 
who studied the regularity of u and of the free boundary d{\Vu\ < n} in the case of a square 
EB, a regular m and an irregular polygon [23] . The key observation is that the problem 0 
is equivalent to the following obstacle problem : 


mm < - 


(5.2) 


| -J \S7u\ 2 dx — J udx : u G Hq(D), u<ndojj j, 

where by doD : D —>► [0, +oo) we denote the distance function 

ddD^x) = min{|x - y\ : yedD}. 

For a general domain D cK d the equivalence of 0 and ( |5.2| ) was proved by Brezis and Sibony 
[3] and is based on a bound on the Lipschitz constant of the solution of (5.2). In fact since we 
have that 

|u G Hq(D) : |V?i| < fcj C 17/ G Hq(D) : u < ftdao j, 


it is sufficient to prove that the solution of the obstacle problem (5.2) is ^-Lipschitz. This may 


not be true in the case of a general force term / when large oscillations can be produced directly 
by / and so the equivalence of the two problems does not hold in general (see m- 
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The regularity of u is due to several authors. In [4j Brezis and Stampacchia showed that 
u G C 1,a (D), for every a G (0,1). The optimal regularity u G C l ^{D) was obtained by Caffarelli- 
Riviere [9], Evans M and Wiegner [24]. In the case of planar domains with boundary which is 
a union of C 3 curves the free boundary was characterized by Caffarelli and Friedman [8]. 

In this case the optimal reinforcement 0 is supported on the set {u = kcLqd} and satisfies 
the Hamilton-Jacobi equation 


V0 • \7dgjj + 0A cIqd H-b A cIqd = 0 on \u = kcIqe)}. 

k, 

This fact was used by Brezis in [2] to obtain 0 explicitly and also to show that 0 is continuous. 
We notice that the regularity of 0 for a general nonconstant / is an open question. We conclude 
this section with an example in which the optimal reinforcement 0 can be easily computed; other 
cases are treated numerically in Section [6j 

Example 5.1 (The radial case). Let D be the unit ball of R d and f — 1. We are then dealing 
with the optimal torsion problem. Passing to polar coordinates gives the optimization problem 

mini- [ r d ~ l \u\ 2 dr + ™ ||?/||^o — [ r d ~ l udr : 7i(l) = oj, 

^ 2 Jo 2 dujd Jo J 


where uo d is the measure of D. After some elementary calculations we find that the solution u 
is given by 


u(r ) 


a 2 — r 2 


2d 


+ 


;(! ~ r ) 


a( 1 — a) 
d 


if r G [0, a] 
if r G [a, 1] 


where a = a m is a suitable number in ]0,1[. Optimizing with respect to a we obtain that a m is 
the unique solution of 

a d+l — {d + l)a(l + —) + d = 0. 

V cu d J 


It remains to compute the optimal function 0. By the Euler-Lagrange equation of 



1 + 0 „d -1 


u 


/12 


dr 


f 

'J dm 


„d— 1 


udr 


we obtain the expression of 0, given by 

0{r) — — -1. 

(dm 


6 Approximation of the free boundary problem 

6.1 A discrete convex constrained formulation 

We introduce in this section, a convergent discretization of the problem 

min {/) (|l Vu l 2 ~ ^“) dx + 5ll V “ll°° : ueH o( D )} ( 6 - 1 ) 

for a domain D cR d . We focus our description on the case of the dimension 2. It is straight¬ 
forward to generalize this approach to the three dimensional case. Our approach relies both on 
classical techniques in approximation of linear partial differential equations and on a suitable 
reformulation of this non-smooth convex energy. Assume r d is a given mesh made of simplexes 
which approximates the computational domain D. In the context of P\ finite elements, we 
consider iL 1 (T)/ l ) the space of functions which are globally continuous and piecewise linear on 
every simplex of the mesh T d . Analogously to the continuous framework, we denote by H^D^) 
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the functions of H l (Dh) which vanish at every boundary point of r^. We define and M\ 
respectively the stiffness and the mass matrices associated to this discretization, more precisely 
we have 


K h 



X7uj • Vun dx 


ho 


M h = 



U{Uj dx 


h3 


where the family (iq) is a basis of the the vector space H^Dh). Thus if we denote by A the 
coordinates of a function of H^Dh) in this basis, we can introduce the discrete E m ( A) energy 
associated to problem (6.1) defined by 


l -\ T K h \_v T f M h \ + ^ sup{(^A)f + (K y X)f} 


where the exponent T denotes the transpose operator, the coordinates of the vector Vf are given 
by the evaluation of the fixed function / on every vertex of the mesh and the matrices K x and K y 
are the linear operator associated to the evaluation of partial derivatives d x and d y in H l (Dh)- 
Notice that K x and K y are not square matrices since their first dimension is indexed by the 
number of simplexes of r^. In order to get rid of the sup term which makes the original energy 
non smooth, we classically introduce a new variable t G M and the new quadratic constraints 

{K x \fi + (K y \)? < t (6.2) 


where l is an index associated to the set of the simplexes of r^. In this context, we define a new 
equivalent smooth convex constrained optimization problem which is the minimization of the 
energy F m (\,t) defined by 


m 


t ) — -A Kh\-Vj Mh A + —t 


for admissible (A,£) which satisfy the constrained (6.2). 


(6.3) 


6.2 Numerical results in 2D and 3D 


In order to study the numerically solution of problem |6.3| under constraints 6.2, we used a 


classical interior point solver provided by the so called IPOPT library available from COIN-OR 
(http://www.coin-or.org). Since the cost function and all the constraints are of quadratic 
type, we provided an exact evaluation of the Hessian matrix with respect to our variable (A, £). 
Moreover, noticed that due to this quadratic formulation we can easily define a fixed pattern 
for the Hessian matrix which allows, through factorization techniques, to speed up dramatically 
the computation of the descent direction at every step. We present below several experiments 
based on this numerical approximation. 

We first focused our study to the example of the unite disk. We discretize this set by a mesh 
made of 2 x 10 4 triangles. This complexity leads to an optimization problem of an approximate 
number of unknowns of the same order and under 2 x 10 4 constraints. We represent on the first 
column of Figure[l]the optimal graphs obtained for m — 0, 0.1, 0.5. We added on the graphs the 
analytically computed radial solution and observe an almost perfect match between our solution 
and the theoretical one. On the second column of Figure [lj we plot the graph of the gradient of 
the solution and for m — 0.1, 0.5. We underlined the free boundary set by adding a bolted line 
along its boundary. We reproduce the same experiments for three other geometries on Figures 
[2j[3]and[4| As expected, we observe in the three cases that for m large enough the free boundary 
of the problem converges to the singular set of the distance function of the domain. 

Finally, to illustrate the efficiency of this approach even in a large scale context, we approx¬ 
imated the optimal profiles on three dimensional problems. In Figures [5] and [6] we discretized 
a cube and a dodecahedron by meshes made of approximately 10 6 simplexes for m = 5. The 
IPOPT library was able to solve such problems in several minutes on a standard computer. We 
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plot in Figures [5] and [ 6 ] different cuts of the graph of the modulus of the gradient of the solution. 
As in the two dimensional case we observe a convergence of the free boundary to a complex 
geometrical set which is very close from the singular set of the distance function. 
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Figure 1: Optimization on a disk for m = 0,0.1, 0.5 
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Figure 2: Optimization on a square for m = 0,0.1,0.5 
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Figure 3: Optimization on an ellipse for m — 0, 0.1, 0.5 
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Figure 4: Optimization on a treffle for m = 0,0.1, 0.5 
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Figure 5: Tomographic view of the optimal solution defined in a cube for m — 5. The free 
boundary is visible in light color and is very close to the cut locus of the cube 
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Figure 6: Tomographic view of the optimal solution defined in a dodecahedron for m — 5. The 
free boundary is visible in light color and is very close to the cut locus of the dodecahedron 
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